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ABSTRACT 

Planets that form around stars born in dense stellar environments such as associations, open clusters, 
and globular clusters are subject to dynamical perturbations from other stars in the system. These 
perturbations will strip outer planets, forming a population of free-floating planets, some of which 
will be tidally captured before they evaporate from the system. For systems with velocity dispersion 
a ^ 1 kms~^, Jupiter-mass planets can be captured into orbits with periods of 0.1 — 0.4 days, which are 
generally stable over ^ 10^ years, assuming quadratic suppression of eddy viscosity in the convective 
zones of the host stars. Under this assumption, and the assumption that most stars form several 
massive planets at separations 5 — 50 AU, I estimate that ~ 0.03% of stars in rich, mature open 
clusters should have extremely close-in tidally captured planets. Approximately 0.005% of field stars 
should also have such planets, which may be found in field searches for transiting planets. Detection 
of a population of tidally-captured planets would indicate that most stars formed in stellar clusters. In 
typical globular clusters, the fraction of stars with tidally-captured planets under these assumptions 
rises to ~ 0.1% - in conflict with the null result of the Hubble Space Telescope transit search in 47 
Tuc. This implies that, if the quadratic prescription for viscosity suppression is correct, planetary 
formation was inhibited in 47 Tuc: on average < 1 planet of mass > Mjup (bound or free-floating) 
formed per cluster star. Less than half of the stars formed solar-system analogs. Brown dwarfs can 
also be captured in tight orbits; the lack of such companions in 47 Tuc in turn implies an upper 
limit on the initial frequency of brown dwarfs in this cluster. However, this upper limit is extremely 
sensitive to the highly uncertain timescale for orbital decay, and varies by four orders of magnitude 
depending on the choice of prescription for the suppression eddy viscosity. Therefore, it is difficult to 
draw robust conclusions about the low-mass end of the mass function in 47 Tuc. 
Subject headings: planetary systems - binaries: close, eclipsing - stars: low-mass, brown dwarfs, 

mass function - globular clusters: general, individual: 47 Tuc - open clusters and 

associations: general 



1. INTRODUCTION 

The discovery of an extrasolar planetary companion 
orbiting a mere 0.05 AU from its host star 51 Peg 
ijMavor fc Ouelod 119951) originally came as a surprise, 
given the seeming impossibility of forming a giant planet 
so near to its parent star. However, the theorists quickly 
recovered, and, drawing on older works which studied 
interactions between protoplanets and th eir natal disks 
^Coldreich TremanufjlTflMlWaTdiriflSa. invoked plan- 
etary migration as the mechanism for delivering the com- 
panion to 51 Peg from its birthplac e to its c urrent ob- 
served position (e.g., iLin. Bo den heimer. fc Ri chardson 

Since this time, considerable progress in the the- 
ory of planetary migration, combined with the detec- 
tion of a large number of additional giant planets on 
short-period orbits, has left little doubt that this pro- 
cess is important in shaping the architecture of plane- 
tary systems. However, many questions remain. For 
example, the mechanisms by which planetary rnigra- 
tion is stopped are uncertain fKuchner fc Lccar' "200^ 
[Lin. Bodenheimer. fc Richardson 1996; Trilling et al. 
11998^ . These mechanisms have the additional burden 
that they should explain both the observed 'pile-up' 
of planetary companions at periods of P = 3 d (i.e. 
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iKuchner fc Lecail 12002*) . an d the n ewest discovery of a 
planet at P = 1.2 d (KonackfeiaOllOOS; Sassclov 2003). 
In addition, it is clear from examples of extrasolar gi- 
ant planets orbiting close to their supposed birthplaces, 
as well as our own solar system, that orbital migration is 
not always so efficient. A satisfactory explanation for the 
causes of these variations in the efficiency of migration is 
lacking. 

It seems plausible that many, if not the majority, of 
disk stars were initially formed in stellar systems - loose 
associations or open clusters that have since dissolved. 
Evidence for thi s comes in part from st udies of the nearby 
moving groups l|de Zeeuw et alJll999|) . Planets orbiting 
stars in stellar systems are subject to numerous effects 
resulting from the dense stellar environment. Planets 
orbiting at larger distances from their parent stars are 
generally more sensitive to these effects. Therefore, the 
structure of planetary systems will be the result of a com- 
plex interplay between the local properties of the system 
that drive planetary formation and migration, and the 
non-local effects that arise from the star's environment. 
Thus, if it is indeed the case that most stars formed in 
stellar clusters, planetary systems cannot be understood 
as isolated systems: a complete picture of planetary for- 
mation and evolution must consider both local and non- 
local effects. 

Several authors have considered the effects 
of dense stellar environments on the formation 
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and surv ival of planetary systems llArmitagd 

IBonnelL .Smith. Davies. Hnrn^ 

Davies fc S igurdss oiJ |200H iHurlev fc Sharal 
Smith fc B onncll 2 001[1 . From these studies it is clear 
that, if distant planets can form in stellar clusters, they 
are likely to be stripped from their parent stars, resulting 
in a free-floating population of planets t hat will slowly 
evapo rate from the system. Recently, IBonnell et al.l 
(|20,0jD suggested that freely-floating brown dwarfs (BD) 
in globular clusters will occasionally be tidally captured, 
resulting in a detectable population of extremely close-in 
BD companions to main-sequence stars if the frequency 
of freely-floating BDs is large. Similar considerations 
can also be applied to planets in stellar systems. 

Here I consider the formation and stability of ex- 
tremely close-in giant planetary companions to stars via 
tidal capture of free-floating planets. In ^\ I summarize 
the dynamical processes of ionization, evaporation, cap- 
ture, and tidal decay, and their relevant timescales. In 
SjSl I give a crude estimate for the frequency of extremely 
close-in giant planets as a function of the parameters of 
the stellar system. In 21 I consider tidally-captured BD 
and planetary companions to stars in globular clusters, 
and interpret the nu ll result of the 47 Tuc transit search 
ijGilliland et al.ll200'0 '). Finally, in I discuss possible 
implications and prospects for the detection of such a 
population of planetary companions. 

2. DYNAMICAL PROCESSES 

It is not clear if planets can form around stars in stel- 
lar systems at all. In particular, in rich clusters with 
N > 10^ members, high-mass stars can generate an ul- 
traviolet radiation field that is sufficiently intense to pho- 
toevaporate protoplanetary discs in a few hundred thou- 
sand years l)Armitagdl2000() . However, planet formation 
might still be possible if there exists a substantial delay 
between high and low-mass star formation. This fact, 
combined with our incomplete understanding of the effi- 
ciency of planetary migration, make it difficult to predict 
a priori the frequency and distribution of planets orbiting 
stars in stellar systems. For the purposes of discussion, 
I will sweep all of these uncertainties under the rug, and 
simply consider the evolution and consequences of a sub- 
stantial population of long-period planets around stars 
in stellar systems. 

Binaries in stellar systems evolve under random en- 
counters with other stars in the syst em. This evo lution 
proceeds according to Heggie's law ()Heggielll97,^ . such 
that soft binaries tend to get disrupted by encounters. 
A binary is soft when its binding energy is equal to the 
mean kinetic energy of the stars in the system. Thus 
there is a critical separation above which binaries will 
tend to be disrupted. 
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where nip is the mass of the secondary, M* is the mass of 
the primary, (M) is the average mass and a is the velocity 
dispersion of the cluster stars, respectively. Planetary 
companions to solar-type stars (wp ~ .^jup and ~ 
Mq) in typical open clusters ((M) ~ 0.5Mq and a ~ 
1.5km s~^ are will tend to get disrupted if a > 0.4AU, or 
P > lOOd. 



The timescale for disruption of a planet is approxi- 
mately the time it takes for the random encounters to 
change the planet 's energy by an a mount equal to its 
binding energy (Bi nnev fc Tremaind HQS?) , 



tdis = (1 + g) 



(2) 



Here q = mp/M,,, v is the stellar number density, and 
InA ~ ln0.4iV is the Coulomb logarithm. The number 
density and velocity dispersion within a stellar system 
can vary by several orders of magnitude. For simplicity, 
here I will adopt these properties averaged over the re- 
gion interior to the half-mass radius of the cluster. 
The mean stellar density in this area is simply 
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For a wide range of equilibrium stellar systems, the half- 
mass radius can be related to the total mass and velocity 
dispersion of the system by ( ^Spitzer 1987.) , 
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I will adopt this relation throughout unless stated oth- 
erwise. For an average stellar density v ^ lO'^ pc~^ 
and velocity dispersion a ~ 1.5 kms^^, the above re- 
lation yields N ~ 7600. For the stellar parameters given 
above, tdis ~ 0.5 Gyr(a/2.5AU)^^. Thus most planets 
with a > 2.5AU in an open cluster will be disrupted over 
the cluster's lifetime (~ 0.5 Gyr). 

Since planets are typically liberated due to the accu- 
mulated effects of distant encounters, free-floating plan- 
ets will initially have velocities similar to the stars, and 
will not escape the system immediately (Hu rley fc Sharal 
[2003- However, equipartition of energy will slowly drive 
the velocity dispersion of free-floating planets to the 
point where many of the planets have sufficient veloc- 
ities to evaporate from the system entirely. Assuming 
that the stellar component of the system dominates the 
mass, the characteris tic timescale for this process is the 
relax ation timescale (jBinnev fc Tremaindll987t ISpitzeJ 

inia), 

^3 
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For the parameters adopted above, this yields tfix ~ 
0.03 Gyr. 

Freely-floating planets will occasionally pass suffi- 
ciently close to another star to raise a significant tide. 
If the energy required to raise this tide is larger than 
the relative energy of the planet and star at infin- 
ity, then the encounter will lead to a bound sys- 
tem (iFabian. Pringle. fc Reegll975l:ILee fc Ostrikeill986l: 
IPress fc Teukolskvlll977ir This leads to the requirement 
that the planet mus t pass within a minimum dis tance 
Ocap to be captured ijFabian. Pringle. fc Reeslll975H . 
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For q = 10-3, - Mq, i?. 

Thus the planet must pass within two stellar radii to be 
captured. 
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Fig. 1. — The lines indicate regions in the plane of planet 
mass and period where various processes arc important, assuming 
solar-mass, solar-radius primaries. Planets can be captured into 
orbits the left of the solid line. Planets captured to the left of the 
long-dashed-dot line result in mergers. Planets to the left of the 
dashed (short-dashed-dot) line are stable over 0.5 Gyr, assuming 
quadratic (linear) suppression of eddy viscosity. Planets to the left 
of the dotted line have radii that exceed their Roche-lobe, assuming 
a planetary radius of i?jup. Planets to the right of the dashed line 
residing in a cluster with velocity dispersion of ~ 1.5km s~^ are 
soft, and tend to be disrupted by perturbations from other cluster 
stars. The lightly-shaded region indicates planet parameters than 
can result in stable captures, assuming quadratic viscosity suppres- 
sion. The solid points show the mass and period known extrasolar 
planetary companions. 



The timescale for tidal capture is, 
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where I have assumed that GM^j^a^ ~ adis/Q ^ Ocap, 
which is vahd for the cases considered here. I find 
^enc — 500 Gyr(acap/2i?0)^^. Assuming every star 
forms Np planets that are liberated from their par- 
ent stars essentially immediately, the fraction of stars 
with tidally captured planet can be crudely estimated as 

Planets passing within acoii — Rp + R* of a star will 
physically collide with it, resulting in a merger. This 
therefore sets an absolute lower limit on the period of a 
tidally captured planet. For gaseous objects supported 
only by degeneracy pressure (including giant planets and 
brown dwarfs), Rp ~ O.li?©, roughly independent of 
mass. For ~ Rq, flcoii — 0.005AU, or a period of 
P ~ 0.13 d. This also sets a lower limit on the mass of 
a planet that can be captured, as very low-mas planets 
must pass within the collision separation in order to raise 
a sufficient tide to be captured. This limit is found by 
equating acoii with Ocap: 



> 



i?*cr / ^ Rp 
GM^ V R* 



(8) 



where I have assumed q ^ 1. For the fiducial parameters, 
this is q > 2 X 10"^ 

Low-mass planets may also be tidally captured suffi- 
ciently close to their parent stars that their radii exceed 
the Roche limit. The limiting separation for Roche lobe 
overflow is. 
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where, for the last relation, I have assumed Rp — O.li?©. 
The Roche limit separation exceeds Ocoii for g > 2 x 10^^^. 
Planets that are captured within their Roche limit will 
lose mass to their parent star. This mass loss is accom- 
panied by a transfer of angular momentum, resulting in 
an outward migration of the planet. Thus a planet ini- 
tially captured inside its Roche limit can still result in 
a stable system. The lower the mass of the planet, the 
farther it must migrate to halt mass transfer. Very low 
mass captured planets may not survive at all. 

Planets captured into very close orbits will quickly 
circu larize, on a timescale of ~ 10^ years l|R.asio et alJ 
119961) . and thereafter be subje ct to tidal d ecay. The 
timescale for tidal decay is (e.g. iRasio et aLl irQQG). 
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where Mcz is the mass of the convective zone; for solar- 
type stars, Mcz ^ Q.02Mo. Here tc is the convec- 
tive eddy turnover timescale, which I will approximate 
as tc \Mr.y^0.2Rl /3Lt: Y/^ , where L^, is the luminos- 
ity of the star (see, e.g.. iRasio et al.lfl996|) . For solar- 
type stars, tc ~ 20 d. Thus for tidally-captured plan- 
ets, P ^ tc, and the largest eddies can no longer con- 
tribute to the total viscosity, resulting is a suppression 
of the viscous dissipation. Therefore / is less than unity. 
Ho wever, the correct form for / remains controversial 
(seelGoodman fc OhI ljl997|) for a discussion). Following 
iSasselovl ()2003D . I will consider both linear suppression 
of the eddy viscosity (|2ahn 1989), 
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as w ell as quadratic suppression ijGoldreich fc KeelevI 
[1973), 

fQ=iT^^ ' ^«*c- (12) 
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Although quadratic suppression is theoretically better 
motivated, linear suppression appears to be in better 
agreement with the obse rved timescale of tidal circular- 
ization in close binaries l|Goodman fc Ohlll997f) . 

For tidally-captured planets with Ocap — 2i?0, the de- 
cay timescale for linear suppression is tdcc — 0.1 Gyr, 
whereas for quadratic suppression, tdec — 80 Gyr. There- 
fore the choice of prescription has an enormous effect on 
the number of surviving tidally-captured planets: one ex- 
pects a negligible number of surviving planets for linear 
suppression, whereas planets are stable over much longer 
than ^ Gyr for quadratic suppression. 

The angular momentum lost from the tidal decay of the 
planet's orbit will be transfered to the star, spinning it 



4 



TIDAL CAPTURE OF PLANETS 



up. Assuming the star is initially slow ly rotating, spin-u p 
of the star will occur on a timescale ijRasio et al.lll99(t) , 



-tdc 



^orb 



q 



tdc 



(13) 



where /* and lorb are the moments of inertia of the star 
and orbit, and fc^ = For a = Ocap ~ 2i?*, and 

fc2 = 0.08, tsu ~ 20(10~^/(/)tdec- Thus companions with 
q < 0.02 will decay before they spin-up the star suffi- 
ciently to synchronize the star's spin period with the or- 
bital period. Conversely, companions with q > 0.02 will 
synchronize before they decay completely. This equilib- 
rium is stable only if the orbital angular momentum of 
the planet is more t han three times the spin angular mo- 
mentum of the star l)Hutl 19801) . This leads to a minimum 
mass ratio for stability. 



q > 3fc2 



(14) 



Thus companions with q < 0.06 cannot reach a stable 
equilibrium, and will always decay. More massive com- 
panions may reach a stable equilibrium. 

The question of orbital stability via synchronization is 
further complicated by the fact that the star will likely be 
spun- up as a result of the capture process itself^, and thus 
may be fairly rapidly rotating. This will not only affect 
the viscous dissipation mechanisms, but will also gener- 
ally shorten the timescale for synchronization. Given the 
uncertainties involved in determining the effects of these 
various process, and since I am primarily interested in 
planetary-mass companions, for which synchronization is 
likely less important, I will furthermore ignore synchro- 
nization entirely, and assume that all orbits decay on the 
timescale given by eq. (|10|) . However, I note that this 
issue has direct bearing on the discussion of the stability 
of tidally-captured BD companions to stars in 47 Tuc in 

SI . 

Figure ^summarizes the regions of parameter space in 
the companion mass-period plane where ionization, tidal 
capture, and tidal decay are important, as discussed in 
this section. For this figure, I have assumed as before 
= Mq, (M) = O.5M0, = Rq^ a = 1.5 kms'^, 
and v = 10'^ pc^^. In addition, I have assumed a fixed 
companion radius of Rp = i?jup, and a system lifetime 
of 0.5 Gyr. Assuming quadratic suppression prescription 
of tidal decay, there is a small sliver of parameter space 
in the period-mass plane of where planets can be tidally 
captured and are stable to tidal decay over typical open 
cluster lifetimes. For linear suppression, there is no re- 
gion of parameter space where tidally-captured planets 
are expected to be stable on such timescales. Also shown 
are known extrasolar planets. Clearly tidally-captured 
planets occupy a region of parameter space that is rela- 
tively disjoint from known companions, and would likely 
to easily identified as the result of capture (rather than 
migration). 

3. FREQUENCY ESTIMATES 

The frequency of tidally-captured planets in a stellar 
system of a given age depends on the competing effects 
of ionization, evaporation, capture, and orbital decay. 

^ I am indebted to I. Bonnell for pointing this out. 



The exact number will depend on the balance between 
the timescales of these various dynamical effects. This 
will in turn depend on the dynamical properties of the 
cluster at all times in its evolution, and thus the pro- 
cesses of relaxation, mass segregation, mass loss, etc., 
should all be considered. Ideally, the most robust way to 
accomplish this is through detailed N-body simulations. 
However, such a study is outside of the scope of this 
paper. Furthermore, such a detailed study is perhaps 
not warranted, as the number of tidally captured planets 
will also depend critically on the initial frequency and 
distribution of planets formed around stars in the clus- 
ter. Given that a general theory of planet formation is 
lacking, this can only be a wild guess at best. I therefore 
simply provide only a crude estimate for the frequency of 
tidally-captured planets. This estimate should be good 
to an order-of-magnitude, and should serve to elucidate 
the dependence of the frequency on the gross parameters 
of the stellar system and the input assumptions. 

I consider only average properties of the stars and plan- 
ets in the stellar system. The dynamical properties of 
the stellar system are specified by the average number 
density v and velocity dispersion a inside the half-mass 
radius, and the average mass of the stars, (M). The 
total number is then set by equations lO and Q). Un- 
less otherwise stated, I assume (A/) = 0.5A/o, which is 
roughly appropriate for mass functions observed in the 
field , young cluster s, open clusters, and globular clusters 
(see lChabriedl2003l) . Here and th roughout, I a dopt radii 
and bolometric luminosities from Allen (1976') and con- 
vective zone masses from Pinsonneault, DcPov. & Coffe^ 
ll20f)l . 

For my fiducial calculations, I assume that every star 
in the system originally has four companions with mass 
Mp = Mjup, with a distribution that is uniform in logo 
between 5 — 50AU. This choice is primarily motivated 
by the distribution of massive planets in our own solar 
system, but is not inconsistent with the distribution of 
extrasolar planets detected via radial velocity surveys. 
These surveys are only just becoming sensitive to plan- 
ets around solar-type stars at 5AU; however extrapola- 
tions based on current samples indicate that the frac- 
tion of stars with Jupiter-mass pl anets at periods longer 
than this may be quite large (e.g.. lLineweaver fc Grethe'3 
I2003|) . For the same distribution of a but other choices 
for the frequency of massive planets, one can simply scale 
all results by {Np/A), where Np is the average number of 
planets per star. 

Free-floating planets in stellar stellar are depleted via 
evaporation and replenished by disruption of bound plan- 
ets"^ . The rate of change in the number Ng of free-floating 
planets per star is thus 



diVfi 
dt 



diV,. 



dt 



dt 



(15) 



Here dN^is/dt is the rate at which planets are ionized. 



diVdi 
dt 



d 
dt 



dK 



dloga 



— 'Pdis(a)dloga 



(16) 



Here dNp/ dloga is the initial distribution of planetary 
companions, which is constant between 5 — 50 AU and 

^ Free-floating planets are also depleted via captures, but the 
rate of captures is generally so small that this can safely be ignored. 
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Fig. 2. — Evolution of the planetary population of stellar clus- 
ters. Every star in the cluster is assumed to initially have A^p = 4 
planets with a uniform logarithmic distribution between 5 and 50 
AU. For other values for the number of planets per star, all curves 
can be scaled as Np/4. Each curve shows the number of planets 
per star as a function of the age of the cluster. The solid curve 
shows the number of planets ionized from their parent stars. The 
dashed curve shows the number of these planets that have escaped 
the cluster; dotted curve shows the remaining free-floating popula- 
tion. The long-dashed curve shows the number of tidally captured 
planets. The short-dashed-dot (long-dashed-dot) curve shows the 
fraction of tidally captured planets that have survived orbital de- 
cay assuming quadratic (linear) suppression of eddy viscosity. The 
vertical line shows trix/?ei the typical timescale for evaporation 
from the system, (a) Assuming a cluster velocity dispersion of 

= 10 kms"'^ and a number density of = 10^ pc 
to a globular cluster, (b) a — ^ ""^ 



appropriate 
3 



priate to a rich open cluster, (c) a = 0.6 kms 
appropriate to a loose association. 



zero otherwise, and Vdis = 1— e" 



-Vtdis(a) 



and u = 10 pc 



is the probabiUty 



that planet with separation a will be disrupted. The rate 
at which planets are evaporated is 



diV,. 



dt 



f ' 



(17) 



where = 0.156 is the evaporation p robability appro- 
priate for test particles l)SDitzeil[l987j) . The number of 
free-floating planets as a function of time can be found 
by integrating eq. (|15|l from t = 0. 

Tidally-captured planets are depleted via orbital decay 
and replenished via new captures of free-floating planets. 
The rate of change in the number Nsn-c of tidally-captured 
planets per star is thus 



dA^s, 



cap 



diVdcc iVff A"s, 



dt 



dt 



dt 



to. 



tdc 



(18) 



I determine the mean capture cross section by averaging 
the width of the region between Ocap and Ocoii, which 
depends on the relative velocity between the planet and 
star, over a Maxwellian (relative) velocity distribution 
with dispersion a. I then evaluate tenc for this cross 
section, and tdoc for the average semi-major axis of the 



captured planets. The number of surviving free-floating 
planets is then found by integrating eq. ((TH|l . 

Note that I have ignored the effects of mass segrega- 
tion. Mass segregation will decrease the rate at which 
planets are captured, as it will partially decouple the dy- 
namical interactions between the stellar and free-floating 
planetary populations. This will generally be a relatively 
small correction, as most planets are tidally captured 
during the first few relaxation times of the cluster. 

Figure 121 shows the resulting evolution of the planetary 
population for three types of stellar systems, character- 
ized by the parameters and a. These are v = 10'* pc~^ 
and cr = 10 kms~^, appropriate to a globular cluster, 
v — 10'^ pc^"^, <T = 1.5 kms~*, appropriate to a rich open 
cluster, and = 10^ pc~'^, cr = 0.6 kms~*, appropriate 
to a loose association. These systems have a total num- 
ber of stars and half-mass radii equal to ~ 7 x 10^, 
rhm = 2 pc (globular cluster), A^ ~ 7600, rhm — 1 pc 
(rich open cluster) and A^ ~ 1500, r/i^ ~ 1.2 pc (as- 
sociation). These values are not meant to be definitive, 
but merely span the interesting range of possible sys- 
tems. Note that r/„„ does not enter into the calcula- 
tions directly, and A^ enters only logarithmically via the 
Coulomb logarithm. 

Due to their high stellar density and velocity disper- 
sions, globular clusters generally have relatively small 
disruption timescales and long relaxation timescales. 
Therefore, planets are liberated from their parent stars 
very quickly, and evaporation is slow, making the dy- 
namics particularly simple: the number density of free- 
floating planets is roughly constant, and d log A"sur/ d log t 
is approximately constant when tidal decay is negligi- 
ble. At the typical globular cluster lifetime of 10 Gyr, 

The evolution of rich open clusters is more complicated 
due to their short relaxation times. From inspection of 
FigureEl it is clear that the majority of planets are tidally 
captured within a time irix/Ce 0.2 Gyr, after which the 
number of surviving planets is approximately constant 
at A^sur — 0.03%. The lifetime of open clusters is uncer- 
tain, but rich clusters have been observed with ages o f 
a Gyr or more ijBurke et al.ll2003at iKalirai et alJ l200lV 
Less massive clusters are probably dissolved on shorter 
timescales. 

Loose associations have even shorter relaxation times, 
but are otherwise similar to open clusters. However one 
important difference is that, because of their smaller ve- 
locity dispersions, planets are captured at larger separa- 
tions. Due to the extremely strong dependence of the 
tidal decay rate on the semi-major axis, this has a pro- 
found effect on the survival of tidally-captured planets 
for the linear prescription for the suppression of eddy 
viscosity: even under such a prescription, one expects 
a significant number of surviving tidally-captured plan- 
ets at a Gyr or less. Loose associations are likely to 
disperse fairly rapidly, on timescales of a few hundred 
million years or less. At 10* years, A'sur — 0.005% 

4. BROWN DWARFS AND PLANETS IN 47 TUG 

Because of their high encounter rates and large num- 
bers of stars, globular clusters are expected to contain 
a large number of systems formed via tidal capture. 
Using the results from the previous section, a system 
with ly ^ 10* pc~"^ and cr 10 kms^* should have 
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~ NNsur = 7 X 10^ X 0.001 ~ 700 systems of tidally- 
captured planets. This process is not limited to planets 
stripped from their parent stars, of course: any plan- 
ets and brown dwarfs formed in isolation will also be 
tidally captured by stars in globular clusters. The ini- 
tial mass function of single objects is observed to be ap- 
proximately flat down to well belo w the hydrogen b urn- 
ing hmit in many you ng clusters ijChabried 120031 but 
see iBricefio et aT]l2CI0l . Such a mass function implies a 
similar number density of stars and brown dwarfs. If this 
mass function is universal then globular clusters should 
also contain a lar ge number of stars w ith tidally-captured 
BD companions ijBonnell et alJl2?)03^ . 

Because of their extremely close orbits, tidally- 
captured planetary and BD companions have a very 
high probability (> 25%) of transiting their parent star, 
making photometric monitoring an id e al wa y of detect- 
ing such companions. IGilliland et alJ l)2000() monitored 
the globular cluster 47 Tuc continuously for ~ 8 d us- 
ing the Hubble Space Telescope, but found no tran- 
sits. This null result can be used to place constraints 
on the initial frequency of planetary companions as well 
as freely-floating planets and BDs formed in isolation 
(jBonncU ct al. 2003). 

Figure El shows the expected number of tidally- 
captured companions per star in 47 Tuc as a func- 
tion of the mass of the companion, for both lin- 
ear and quadratic prescriptions for viscosity suppres- 
sion. I have adopted parameters appropriate to 
47 Tuc, namely i/ = 10* pc'^, cr = 12 kms-\ 
rhm. = 3.9pc, (M) = 0.56Mf7^ and an age of 13 Gyr 



Quadratic 



Linear 



6.5 



jlpiorgovskf "1993; Gcbhardt. Pr yor. Williams, fc Hessej 
ll995iL iParesce fc De Marchi .20001) . Note that here I do 
not assume eq. (@J. FigureOI assumes an initial frequency 
of one object (planet or brown dwarf) per star. Due to 
the extremely high disruption rate, these results are es- 
sentially independent of whether the objects are initially 
free-floating, or bound on orbits with a ^ 1 AU. For the 
quadratic prescription for viscosity suppression, the ex- 
pected frequency of tidally-captured BD (O.OSMq) com- 
panions is A^sur ~ 0.002 companions per star, whereas 
the frequency of tidally-captured Jupiter-mass planets is 
AT^^- 0.0002. 



Gillilan d et all l|2000l) claim that, if / = 0.8 - 1.0% 
of the 34,091 main-sequence stars they monitored had 
1.3-Rjup companions at periods of 3.5d, then they should 
have seen 17 planets. This implies a 95% confidence level 
(c.l.) upper limit of 3//17 = 0.14-0.18% to the fraction 
of main-sequence stars with such companions. This null 
result can also combined with the prediction for the num- 
ber of tidally-captured companions Ngur to determine an 
upper limit to the initial frequency Np of objects per star 

3 / 7'T(3.5d) ^^g^ 



Np < 



17A^surPT(Pcap)' 



Here Vt = {Rp + R*)/a is the transit probability. Figure 
13 shows the resulting 95% c.l. upper limit to the initial 
frequency of objects as a function of the mass of the 
object. The shaded region comes from varying / in the 
range 0.8% < / < 1.0%. 

The null result implies that, on average, each star in 47 
Tuc originally formed < 1 planet with mass Mp > Mjup, 
or < 0.4 planets with Mp > 5Afjup. Less than 50% of 




10-* 10-3 0.01 0.1 10-* 10-3 0.01 0.1 
Mass (Mg) Mass (Mg) 



Fig. 3. — (a) Number of tidally-captured companions per star 
in 47 Tuc as a function of the mass of the companion in solar 
masses, assuming an initial frequency of one per star. A stellar 
density of = 10 pc"'^, velocity dispersion of cr = 12 Itms^^, half- 
mass radius of r;,^ = 3.9pc, and an age of 13 Gyr were assumed. 
The solid curve is for the nominal quadratic assumption of /g = 
(P/27ric)^ for the suppression fg of eddy viscosity, whereas the 
dotted line is for the assumption /q = {P/2tc)'^ adopted by some 
authors. Here P is the period of the planet and tc is the convective 
eddy turnover time. The dashed line is the result for no tidal 
decay, (b) 95% confidence level upper limits to the initial number 
of objects (brown dwarfs or planets) per star implied by the null 
result of the transit search in 47 Tuc (Gilliland ct al. 2000), as a 
function of the mass of the object. Upper and lower curves bracket 
the range of uncertainty in the limits. See text. The lower, middle, 
and upper curves correspond to the limits from the upper, middle 
and lower curves in panel (a), respectively, (c) Same as panel (a), 
but assuming linear suppression of eddy viscosity = {P/2tc). 
(d) Same as panel (b), but for /i. 



stars in 47 Tuc could have formed solar-system analogs. 
While this limit is considerably weaker than the original 
limit placed bv IGillil and et all l|2000(l of ~ 0.1%, it is 
important to emphasize that this original limit applies 
only to planets in close (~ 3 d) orbits. Since planets 
in such small orbits were not formed in situ, and must 
be the result of migration, drawing conclusions about 
the formation of planets at larger separations based on 
the lack of planets at small separations is subject to un- 
certainties in the timescale for orbital migration. If the 
migration timescale were greater than or of order the dis- 
ruption timescale, the planets would have been stripped 
from their parent stars before they had a chance to mi- 
grate sufficiently close to be detectable via transits. In 
contrast, the limit on the frequency of planets implied 
by the lack of tidally-captured planets, although weaker, 
is essentially independent of where they were initially 
formed. In fact, one can conclude the initial number of 
free-floating planets (planets formed in isolation via gas 
fragmentation and collapse) is also less than or equal to 
the number of stars. 
Under a similar set of assumptions as adopted here, 
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but using a more sophisticated s imula tion and neglect- 
ing orbital decay, iBonnell et al.1 l)2003D found an upper 
limit to the BD frequency in 47 Tuc of ~ 15%. Figure O 
shows the upper limit derived using the formalism in ^ 
but neglec ting orbital de c ay. I find 12%, in good agree- 
ment with IBonnell et all l|2flfl,1l) . This suggests that, in 
the absence of orbital decay, the predicted frequency of 
tidally-captured BD companions should be reasonably 
robust. 

For quadratic suppression of eddy viscosity, the 
timescale for orbital decay of tidally-captured BD [q ~ 
0.1) companions is idoc ~ 9 Gyr. Therefore, approxi- 
mately exp(— 9/13) ~ 50% of the captured BD will still 
be present today. Thus tidal decay is significant at the 
current age of 47 Tuc for tidally-captured BD compan- 
ions. Figure 13 shows the upper limit including orbital 
decay. The li mit is revised upward by ^ 2 relative to 
the estimate of IBonnell et al.l l)2003j) . The lack of tran- 
siting BD companions in 47 Tuc implies that the initial 
frequency of BDs relative to stars is < 25%. 

The upper limit of ~ 25% on the relative frequency 
of BDs in 47 Tuc is inconsistent with the observed mass 
functions of most young clusters and th e field, for whic h 
BDs and stars exist in equal numbers l)Chabriedl2003j) . 
although it is consistent with the frequency of BDs in the 
Taurus star- forming region (Briccno et al. 2002). At face 
value, this is an indication that the initial m ass functio n 
compact objects is not universal ( Bonnell et al.l 1200^ . 
Unfortunately, this conclusion relies heavily on the as- 
sumed value for the tidal decay timescale. Different as- 
sumptions for the mechanisms for tidal dissipation will 
give rise to radically different conclusions. For example, 
in Figure |31 1 show the predicted frequency of tidally- 
captured companions under the assumption of linear sup- 
pression of eddy viscosity. The frequencies are always 
< 3x 10~^, i.e. ^ 5 in the entire cluster. This low number 
results from the fact that the tidal decay timescale under 
linear suppression is smaller by a factor of -\/7q ~ 10"''. 
Thus tdec — 4 X 10^ yrs, and the frequency of tidally- 
captured objects is only ~ idoc/^onc 10~^. Obviously 
the resulting upper limits are not very interesting. Even 
if one accepts that the quadratic suppression of eddy 
viscosity is correct, there are still some ambiguities. For 
exampl e, there exist discrepan cies in the quoted form 
for /q: IGoodman fcTThI ifTflflTl) a,nd ISa,sse]ovl ll2nnl use 
the form adopted here, namely fq = (P/27rtc)^, but 
iRasio et"aLl l|1996D adopt /g = /qtt^, i.e. larger by al- 
most an order of magnitude. This factor translates di- 
rectly to the implied upper limit on the relative fre- 
quency, weakening the constraint considerably (see Fig. 
O. Note that this would also weaken the constraint on 
planets in 47 Tuc, but by a smaller factor than the con- 
straint on BDs. 

Thus, given the current uncertainties in the physics of 
tidal dissipation, I conclude that no robust inferences can 
made about the frequency of brown dwarf companions in 
47 Tuc, and any conclusions regarding the universality of 
the initial mass function are probably premature. 

5. IMPLICATIONS AND PROSPECTS FOR DETECTION 

It is clear that, under certain sets of assumptions, ex- 
tremely close-in tidally-captured massive planets should 
exist around stars formed in dense stellar systems. What 



are the prospects for the detection of such companions, 
and what would be the implications of any such detec- 
tions? 

Given their extremely close orbits, the most promising 
method for detecting companions is transits. For a planet 
captured at an orbit of a < 2i?Q, the transit probability 
is Vt ^ 50%. The duration of the transit is 



tT 



P 



(20) 



where i is the inclination. Thus for a central transit, 
the duty cycle (fraction of time in transit) is quite large: 
tr/P ^ 15%. For a solar-type primary, the transits 
will last ^ 50 min and recur every ~ 5 hours. The 
transit depth will be ~ 1%. Given the extremely close 
proximity of the planet to the star, one might worry 
that the planet will induce additional photometric vari- 
ations on the star. The reflected light from the com- 
panion will induce an fractional flux variation of order 
AF/F ^ p{Rp/aY ~ 2 X 10"'^, where p is the geomet- 
ric albedo, and I have adopted p = 2/3 as appropriate 
for Lambert sphere scattering. Although five times lower 
than the signal from the transit, a signal of this magni- 
tude may be large enough to be detected. In fact, de- 
pending on the selection criterion used to choose targets 
for the transit search in 47 Tuc, this may be cause for 
concern in interpreting this result. The planet will also 
induce photometric variations of the star from tidal and 
Doppler effects, b ut these will likely b e small for plane- 
tary companions IILoeb fc Gaudill2n0l . They may, how- 
ever, be significant for more massive companions. Thus, 
barring an unforeseen effects that induce additional pho- 
tometric variability, it should be possible to cleanly de- 
tect tidally-captured planets via transits. 

Because of their short periods and relatively large duty 
cycles, the observational requirements of a transit search 
for tidally-captured planets are generally much less se- 
vere than transit searches for more distant planets. In 
particular, complete coverage of several phases can be 
achieved in less than a few days. This is especially im- 
portant for ground-based transit searches, as the coher- 
ence time of weather patterns is typically several days. 
This proves to be devastating when searching for plan- 
ets with periods of order the coherence time, but nearly 
optimal for very short periods. Furthermore, aliasing is 
minimized because the period of the planet is typically 
less than the duration of a typical clear observing night. 

Tidally-captured planets will also induce radial ve- 
locity variations in the host star, of amplitude K ~ 
340 m s"^(5 hr/P)(Mp/Mjup). For reasonably bright 
sources, this is well within reach of current instrumenta- 
tion. However, the presence of such a close planet may 
induce additional radial velocity variability, making de- 
tection of the planetary signal difficult. Furthermore, 
the extremely low expected frequency of tidally-captured 
planets makes this an extremely inefficient method of de- 
tection, unless multi-object spectrographs capable of pre- 
cision radial velocity measurements become available, so 
that many stars can be monitored simultaneously. 

The systems that offer the best chance of observing 
tidally-captured planets are globular clusters, because 
they have the highest frequencies of captured planets 
(~ 0.025 — 0.1%), and the largest total number of sys- 
tems. Furthermore, due to the large velocity dispersions. 
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planets are captured into very close orbits. This increases 
the a priori probability that the planet will transit its 
parent stars to nearly unity. In ^21 I considered in detail 
frequency and detectability of tidally-captured planets 
in 47 Tuc. Although no transiting planets were found in 
this cluster, it may still be worthwhile to search in other 
clusters. An observing campaign with the Hubble Space 
Telescope need not be as ambitious as that toward 47 
Tuc. 

The expected number of transiting, tidally-captured 
planets in a rich open cluster is only of order unity. 
Therefore, the prospects of detecting s uch planets 
in transit searches towa rd open clusters ijBurke et al.l 
EpOSb; Str eet et al.ll20nl are poor, unless many massive 
planets are formed per star. 

Loose clusters are subject to fairly rapid disruption, 
after which their stars disperse into the general disk pop- 
ulation. It is possible that a large fraction of stars in the 
disk were initially born in such environments. Assum- 
ing a typical lifetime of a loose cluster of ^ 300 Myr, 
then ~ 0.005% of stars in the field should have a tidally- 
captured planet, if all stars in the field were born in as- 
sociations, and planet formation were ubiquitous. Ac- 
counting for a binary fraction of 50%, and assuming 
a transit probability of ~ 40% for planets captured in 
loose clusters, I estimate that 100, 000 stars need to 
be monitored to have a chance of detecting even one 
transiting tidally-captured planet. OGLE monitored 5 
million stars toward the Galactic bulge ov er 45 nights to 
searc h for transiting planets to disk stars (lUdals ki et al.l 
120021) . They only analyze a subset of ~ 52,000 of these 
stars with photometry better than 1.5%. Similarly, the 
EXPLORE project has monitore d ^ 350,000 stars in 
the G alactic plane for 11 nights ijMallen-Ornelas et al.l 
l2003fl . They achieved better than 1% photometry for 
a subset of 37,000 stars, which they visually inspected 
for transits. In both cases, an insufficient number of 
stars were searched to place interesting constraints on 
the fraction of tidally-captured planets. However, be- 
cause of the large number of expected transits, it should 
be possible to detect transiting planets even when the 
scatter is larger than the expected signal. Therefore, 
the optimal (but time-consuming) approach would be to 
search all light curves for transiting planets. For such 
a photon-limited survey, the detection probability scales 
as (PcDDcr, Gould & DcPov 2003), highly favor- 

ing extremely close-in tidally-captured planets. 

It is important to reiterate that all of these frequency 
estimates rely on the assumption that the tidal decay 
timescale is longer than ~ 1 Gyr, and thus that the 
quadratic prescription for the suppression of eddy vis- 
cosity is approximately correct. The only exception is 
planets captured in loose associations. Such planets are 
captured sufficiently far from their parent stars that they 
are stable over ~ Gyr timescales for both linear and 
quadratic prescriptions. However, for ages much larger 
than 1 Gyr, the orbits will rapidly decay under the linear 
prescription (see Fig. 0) . The detection of a planet on 
an extremely close orbit would therefore have important 
implications for the tidal dissipation theory. 

How can tidally-captured planets be distinguished 
from planets that migrated to such close orbits? The pe- 
riod distribution of radial-velocity planets shows a pile- 



up at P 3 d, with a significant lack of planets with 
smaller periods. Tidally-captured pla nets would be eas- 
ily identified in such a distribution. iKuchner &: Lecail 
( 2002) suggest that this distribution arises from the fact 
that the centers of protoplanetary disks are evacuated in- 
terior to separations corresponding to periods of 6 d, and 
that planetary migration thus halts when the planet's 
outer Lindlbad resonance reaches the inner disk edge. 
Although this seems to provide a compelling explana- 
tion for the observed period distribution of radial veloc- 
ity planets, the recent detection of a transiting planet 
at P = 1.2 days ijKonacki et al.ll2003l| generally sug gests 
that this pictur e cannot be universal || Sa,sse]ovl2003f) . In- 
deed, iTrilling e t alJ l)1998|) suggest another mechanism 
for halting migration: planets migrate until they reach 
their Roche limit, at which point they lose enough mass 
to the parent star to balance inward migration. If disk 
dispersal occurs in a sufficiently short time, the planets 
will be left intact in short-period orbits. Such planets 
would generally have similar periods as tidally-captured 
planets. Therefore, additional diagnostics are needed. 

One possible method of distinguishing tidally-captured 
planets from migrating planets is a radius measurement. 
Both known transiting planets have radii significantly 
larger than that of Jupiter. It has been suggested that 
these large radii are due to the fact that these planets 
migrated on a sufficiently short timescale that their grav- 
itational contraction wa s retarded, leaving th em in a per- 
manently inflated state ijBurrows et al.l20()(H) . If this pro- 
cess is universal, then we can expect migrated planets to 
have significantly larger radii than tidally-captured plan- 
ets, which could not have migrated far from their birth- 
places. This assumes that there are no processes that can 
inflate planets once they have been tidally captured. In 
principle, variable tidal forces from the star on a planet 
in an eccentric orbit can deposit enough energy to inflate 
the planet IjGu. Lin, fc Bodenheimeill2003l) . but to have 
a significant effect, there must be some external force 
continuously pumping the eccentricity, since the circu- 
larization timescale for extremely close-in planets is very 
short. 

6. SUMMARY 

Planets which form around stars born in dense stel- 
lar environments are subject to dynamical perturbations 
from other stars in the system. These perturbations will 
generally serve to strip outer planets from their parent 
stars, leading to a significant population of freely-floating 
planets in the system. Further interactions with stars 
in the system drive this planetary population toward 
equipartition, thus raising the velocity dispersion of the 
planets, ultimately resulting in escape from the system. 
However, some planets will survive long enough to be 
pass sufficiently close to another star to be tidally cap- 
tured. 

The stability of tidally-captured planets against orbital 
decay depends critically on the physics of viscous dissi- 
pation of the tide in the convective envelope of the star. 
This processes is extremely uncertain, and different pre- 
scriptions lead to decay timescales that differ by four 
orders of magnitude. For quadratic suppression of eddy 
viscosity, planets on tidally-captured orbits will gener- 
ally be stable for > Gyr, whereas for linear suppression, 
planets will decay very quickly. 
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The frequency of tidally-captured planets depends on 
the competing effects of dynamical stripping, evapora- 
tion, capture, and orbital decay. The timescales for these 
effects depend in turn on the properties of the stellar 
systems, particularly the velocity dispersion cr and the 
number density v. Under the assumption that all stars 
form four Jupiter-mass planets with a uniform logarith- 
mic distribution in semi-major axis between 5 — 50AU, 
and assuming quadratic suppression of eddy viscosity, I 
have estimated the frequency of tidally-captured plan- 
ets in a stellar system as a function of the age, velocity 
dispersion, and number density of the system. 

For loose associations with a ^ 0.6 kms~^ and v ~ 
10^ pc~'^, the frequency of tidally-captured planets after 
a typical cluster lifetime of 10'* years is ^ 0.005%. If most 
stars formed in such systems, this is roughly the expected 
frequency of extremely close-in planetary companions to 
stars in the field. These planets may be found by deep, 
wide- angle searches for transiting planets aroun d disk 
stars (|Mallen-Ornelas et all l20M lUdalski et alJl2002D . 
Detection of such a population would imply both that 
most stars were formed in dense stellar environments, 
and that the timescale for orbital decay is > Gyr. 

For rich open clusters [a ^ 1.5 kms~^ and v ^ 
10^ pc~^), the frequency of tidally-captured planets is 
^ 0.03%. Unfortunately, the prospects for detecting such 
a population of planets in a given cluster are poor, as < 1 
transiting planet per cluster is expected. 

Due to the long relaxation times and high encounter 
rates, the frequency of tidally-captured planets in a typ- 
ical globular cluster (tr ^ 10 kms~^ and v ~ lO'' pc~^) 
is ~ 0.1%. Therefore, one expects several hundred tran- 
siting planets in each globular cluster. 

The relatively high frequency of tidal capture in glob- 
ular clusters makes them excellent targets to search for 
such systems. The null result of the search for transiting 
planets in 47 Tuc can be used to place interesting con- 
straints on the initial frequency of planets and brown- 
dwarfs (either initially free-floating or bound) relative to 
stars. I find that < 1 planet with mass > Mjup origi- 



nally formed per star in 47 Tuc. Less than 50% of stars 
formed solar-system analogs. The initial frequency of 
brown dwarfs relative to stars is 25%. All of these con- 
clusions are under the assumption of quadratic suppres- 
sion of eddy viscosity. If viscosity suppression is much 
less efficient, no interesting constraints can be placed. 

Tidally-captured planets will have orbital periods of 
^ 0.1 — 0.4 d, and can be distinguished from migrat- 
ing planets from their distribution of orbital periods pro- 
vided the mechanism that is observed to halt planetary 
migration at P ~ 3 d is robust. However, if the recent 
discovery of a planet at P ^ 1 d is indeed indicating 
that there exist multiple mechanisms for halting migra- 
tion, then one cannot rule out the hypothesis that an ob- 
served extremely close-in planet was the result of orbital 
migration. However, the radii of tidally-captured planets 
are expected to be similar to their radius at their sites 
of formation, and thus considerably smaller than planets 
that migrated quickly to close orbits. Thus a radius mea- 
surement might provide a diagnostic for tidally-captured 
planets. 

If the majority of stars is formed in dense stellar en- 
vironments, then most planetary systems will be subject 
to a variety of dynamical effects that can alter planetary 
formation, migration, and survival. Therefore, planetary 
systems cannot be understood as isolated systems, and 
the interpretation of the observed properties of extraso- 
lar planets will require the consideration of these effects. 
This study highlights one example of how such considera- 
tions can lead to new observable diagnostics of planetary 
formation. 
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